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I. INTRODUCTION 

The unified model of dark energy and dark matter 
[11, 0, S] is an interesting alternative to account the ac- 
celerated expansion of the universe. It is shown to be 
consistent with first year WMAP and other CMBR data 
0, S Hi I gravitational lensing , supernova data [l, Q , 
cosmic topology [ifl] , gama ray burs ts llll a nd structure 
formation (see discussions in Refs. (l2l. Il3||). The com- 
patibility of the GCG scenario with the neutrino back- 
ground radiation has been studied in Refs. [l^. [Tsj. 

The GCG behavior is described by an exotic perfect 
fluid with the following equation of state P, [2] 
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where A is a positive constant and a is a constant in 
the range < a < 1. The covariant conservation of 
the energy-momentum tensor for an homogeneous and 
isotropic spacetime implies, in terms of the scale factor, 
a, that: 
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where B is an integration constant. That is, at early 
times the GCG behaves like non-relativistic matter, while 
at late times it mimics a cosmological constant. This be- 
havior is maintained even after considering energy den- 
sity perturbations 0, 0] ■ 

In this work, we update the analysis carried out in 
Refs. [3, 0] constraining the parameters of the GCG 
model. We use the bounds from the positions of the 
peaks and troughs of the CMBR power spectrum that 
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have an increased precision from the WMAP Three Year 
Observations and we add the bounds arising from the 
baryon acoustic oscillations (BAO). In this scenario we 
consider the GCG to be a unified dark matter/dark en- 
ergy fluid with no additional cold dark matter. 

We also perform a separate analysis using the full 
CMBR spectrum. In this case, we include an additional 
cold dark matter fluid, and allow for the full perturba- 
tion of the GCG fluid, following Refs. @, 13. We use 
a modified CAMB code to compute the theoretical 
CMBR spectrum and the recently released WMAP5 Five 
Year Observations [l8^ to obtain the likelihoods. We then 
perform a Markov chain Monte Carlo analysis with the 
COSMOMC code [l^, using a MetropoHs algorithm, to 
constrain the GCG parameters. 

In section |ll] we describe the theoretical framework of 
our study; our results are presented in section IIIII Fi- 
nally, our conclusions are discussed in section [TVl 



II. 



CMBR CONSTRAINTS ON THE GCG 
MODEL 



A. Constraints from peak and trough locations 

The CMBR peaks arise from acoustic oscillations of 
the primeval plasma just before the universe becomes 
transparent. The scale of these oscillations is set by the 
acoustic scale, Z^, which, for a flat universe, is given by 
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where r = rj^^ J a~^dt is the conformal time, tq and 
Tis being its value at present and at last scattering re- 
spectively, while cs is the average sound speed before 
decoupling: 



= T;/ / CsdT , (4) 
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where 
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with Pbl p-y being the ratio of baryon to photon energy 
density. 

In an ideahzed model of the primeval plasma, there 
is a straightforward relationship between the location of 
the m-th peak and the acoustic scale, namely Im ~ mlA- 
However, the peak positions are shifted by several effects 
which can be estimated by parametrizing the location of 
the m-th peak, as in Refs. 20, 21, 22] 
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where ^ = is the overall peak shift and SLpm = '^m — ^ 
is the relative shift of the m-th peak relative to the first 
one. Eq. ^ can also be used for the position of troughs 
by setting m = 3/2 for the first one. We use the fitting 
formulae of Refs. [13, HH [H, [H for the shifts and the 
theoretical estimate of the peak locations as described in 
the Appendix. 

The energy density, Eq. can be rewritten as 
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through the definitions As = ^/ pjio P^hQ = [A 
_B)T+° . It follows for the expansion rate 
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where pho and pro are the baryon and radiation energy 
densities at present. We do not include cold dark matter 
in this model, its role being taken by the GCG. 

It is worth mentioning that < As < 1 and that, for 
As — 0, the GCG behaves as dust while for As = 1 it 
behaves like a cosmological constant. The GCG model 
matches the AC DAI model for a = and As = 1. 

Using the fact that Pro/ Pcho = f^ro/(l ^^rO~^bo) a-nd 

Pbo/PchO — f^bo/(l — ^rO — ^bo), 



i/2 = nMH^a-^X^a) 
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Moreover, since — a *(^)^>we get 
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so that 
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where ais is the scale factor at last scattering, for which 
we use the fitting formulae 0, [HI, [H, [ll] : 



where 
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and LOb.rn = ^b,mh^- 

For the position of the peaks we consider the three 
year WMAP measurements, which show a considerable 
precision on the locations of the first two peaks and the 
first trough, namely [2^ : 



Ipi = 220.8 ±0.7 



lr,-2 = 530.9 ±3.8 
412.4± 1.9 
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at la. 

We also consider the bound from the baryon acoustic 
peak [25j . Its position is related to the quantity 
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which takes the value 
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An = 0.469 ±0.017 . 



(16) 
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where zirg = 0.35 is the redshift from the Sloan Digital 
Sky Survey for luminous red galaxies. 

Combining Eqs. (|6]) and l|12p and the fitting formu- 
lae of Refs [iO, [U, [52, [11] , shown in the Appendix, we 
search for the combination of GCG model parameters 
that is consistent with the observational bounds, taking 
also into account the theoretical error estimates given in 
the Appendix. For this analysis we fixed the value of the 
spectral index Ug = 0.963 to the WMAP mean, and we 
used LOro = 4.1532 x 10~^ and Ubo = 0.027 for the present 
energy densities of radiation and baryons respectively. 

Our results are shown in section Hill 



B. Likelihood analysis 

Besides the peaks and trough analysis described above 
we also consider the full temperature and polarization 
spectrum to constrain the GCG model with the latest 
CMBR data from WMAP5. We relax our value for the 
cold dark matter, so that we are now really constrain- 
ing a mixed model, where the GCG will be mainly re- 
sponsible for the dark energy. In order to do this, we 
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obtain the theoretical CMBR spectrum using a modi- 
fied CAMB [13] code to include the GCG evolution in 
the background, Eqs. (l7])-([8]), as well as its perturba- 
tions Using derivatives with respect to conformal 
time and defining the conformal time Hubble parameter, 
TC = a' /a, (where the primes denote differentiation with 
respect to the conformal time) we have 
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where 5 and 6 are the density and velocity divergence 
perturbations for the GCG fluid, and 
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are, respectively, the GCG fluid equation of state and 
speed of sound. We assume that the shear of the GCG 
fluid vanishes. 

We included these equations into the CAMB code [l3| , 
and used the latest WMAP five year data likelihood code 
[l^ to compute the likelihoods for this model. We used 
the COSMOMC ^ code to perform a Markov chain 
Monte Carlo analysis on our parameter space. For sim- 
plicity we flxed the inflationary input (namely, the spec- 
tral index = 0.963), and allowed changes in the pa- 
rameters a and As for the GCG fluid, ujh and lVc for the 
present baryon and cold dark matter energy density, r the 
optical depth and h the Hubble constant. We flxed the 
radiation value at the present to be curO = 4.1534 x lO"'"^. 

We present our results in the next section. 
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FIG. 1: Contour plots of the first two Doppler peaks and first 
trough locations, as well as for the baryon acoustic peak, in 
the {Q,rn,h) plane for the GCG model, with ns = 0.963 for a = 
(top graphic) and a = 0.2 (bottom figure). Full, dashed, 
dot-dashed and dotted contours correspond to observational 
bounds on Ipi, lp2, Idi (Eqs. (llSp ) and the baryon acoustic 
peak (Eq. (|17|) ). respectively. The box indicates the bounds 
on h and flmh^ from the WMAP five year data. Darker and 
lighter shaded areas correspond to the allowed regions with 
and without BAO, respectively. 



III. RESULTS AND DISCUSSION 

A. Peak and trough constraints 

Fig. [T] shows the contour plots of the two flrst Doppler 
peaks and the flrst trough locations, as well as the baryon 
acoustic oscillations (BAO) peak, in the {flm, h) plane for 
the GCG model. Since we do not have cold dark matter 
in this model, the Qm parameter is really an estimate 
of the "matter" component of the GCG with the baryon 
energy density, namely through the relation to As given 
by 

il,n = l-^r- As{l - fib - ftr) . (22) 

We see that the dependance of the BAO contour on 
h is very slight, so the BAO constraint works as a lower 
bound on the value of flm- The boxes depicted in the 
figures are the 1-sigma bounds on h and flmh^ from the 
WMAP five year data fit to a ACDM model, as indicated 



in the caption. We can see that for a — (top figure) the 
results are compatible with the observations, as expected; 
notice also that the second peak is not measured with a 
sufficient accuracy to further constrain the results from 
the other bounds. 

For a = 0.2 (bottom figure), we see that the BAO 
bound combined with the first peak bound pushes h to 
smaller values, making the model less viable. Higher val- 
ues of a are almost ruled out from the observations (see 
below). 

Fig. [5] shows the same bounds in the {As, a) plane. 
We now fix /i = 0.719, still using Ug — 0.963. We see 
that the first peak and trough constraints allow large As, 
As — 0.9, and a < 0.4 (again, the second peak location 
does not noticeably change our results). Introducing the 
BAO constraint, however, severily affects these results, 
requiring a < 0.1 and As ~ 0.77. 

In comparing these results with , we have to consider 
that we are now including the estimated theoretical errors 
for the fitting formulae |23|. Taking this into account 
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FIG. 2: Contour plots of the first two Doppler peaks and 
first trough locations, as well as for the baryon acoustic peak, 
in the {As, a) plane for the GCG model, with ns ~ 0.963 
and h = 0.719. Full, dashed, dot-dashed and dotted contours 
correspond to observational bounds on, respectively, Ipi, lp2, 
Idi (Eqs. (HU) and the baryon acoustic peak (Eq. (|17|l '). 
Darker and lighter shaded areas correspond to the allowed 
regions with and without BAO, respectively. 
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FIG. 3: The spectrum for the GCG model, compared with 
a AC DM model (dashed curve), for a = 0.5, ns = 0.963, 
T = 0.087 and h = 0.719, for As = 0.7, 0.8, 0.9, 0.99, 0.999, 
from top to bottom, respectively (at the plateau). 
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shows that the new WMAP results make a considerable 
improvement in constraining the parameter space of the 
model. The BAO bounds add a significant contribution 
to these constraints, further reducing the allowed region 
of parameter values. 
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B. Likelihood analysis 

We start recalling that in this approach, we use the 
full set of perturbation equations for the GCG, so we 
are mainly considering the GCG fluid to act as the dark 
energy component. That is, we allow a separate cold dark 
matter energy density as a free parameter in our model. 

In Figs. [3]and[3]we show the WMAP normalized power 
spectrum computed with CAMB . We used the WMAP 
five year mean values for the spectral index, ns = 0.963, 
the optical depth to the last scattering surface, r — 0.087, 
and the Hubble constant, h — 0.719. We can see that 
the dependance of the GCG model on As is stronger 
than on a. The first and second peak move to the left 
as we decrease As or as we increase a, in consistency 
with the results from the previous sub-section. For the 
same parameter changes, the third peak moves to the 
right. The ratio between the first peak amplitude and 
the plateau gets smaller with decreasing As or increasing 
a, due to the integrated Sachs- Wolfe effect (ISW) [l,[ll]. 
The first and second peak intensities increase for greater 
values of As and smaller values of a, while the plateau 
and third peak decrease. We can see that the model 
favours As values close to one and small a values. 

Fig. [5] shows the marginalized posterior probability 
densities for a and As- We find the bounds 

a < 0.25 , 
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FIG. 4: The spectrum for the GCG model, compared with 
a AC DM model (dashed curve), for As = 0.99, ns = 0.963, 
T = 0.087 and h = 0.719, for q = 1, 0.5, 0.2, 0.1, 0, from 
top to bottom, respectively (at the plateau). 

0.93 < As <l , (23) 

at 2a. 

The obtained mean values for the other parameters are 
LUc = 0.0935, ojb = 0.0225, r = 0.70 and h = 0.74. 

The difference in the bounds of As obtained by the 
previous approach and this one is easily understood. In 
the first case we do not consider the dark matter com- 
ponent, so the GCG is modeling both, dark matter and 
dark energy. In the second approach, dark matter is in- 
cluded, so the GCG contribution to this component is 
not the dominating one. This tends to produce a higher 
lower bound for As- 

IV. CONCLUSIONS 

In the present work we have considered the position of 
the first two peaks and trough from WMAP three year 
data and other CMBR experiments that best fit WMAP 
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are fitting coefficients, ujb = il^^^ 



and 



[Zls) I Pm{zis) 



(A.3) 
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is the ratio of radiation to matter at decoupling. A suit- 
able fitting formula for zis is .27.] 



FIG. 5: Marginalized posterior probability densities for the 
model parameters. Dotted lines are the likelihood functions. 
The vertical dashed lines mark the la and 2a confidence lev- 
els. 



zis = 1048[1 + 0.00124wj-°-^3*^][l 



where 



(A.5) 



five year data for h and ns- Following the analysis of 
Refs. the new WMAP CMBR data shows that the 

GCG model is compatible with observations for a < 0.4. 
The inclusion of the BAO bound further constraints the 
model, the allowed region of the GCG parameters being 
a < 0.1 and As ~ 0.77. 

Secondly, we considered the full CMBR spectrum al- 
lowing the full perturbation of the GCG fluid. We 
used a modified CAMB code to compute the theoret- 
ical spectrum and the COSMOMC code to perform a 
Markov chain Monte Carlo analysis with a Metropolis 
algorithm. The bound on the GCG parameters obtained 
were a < 0.25 and As > 0.93, at 2a. 
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APPENDIX 

We have used the analytic approximations for the 
phase shifts found in Refs. [l^H^I, which we reproduce 
here for convenience. The overall phase shift is given by 



<f = (1.466 - 0.466ns) 
where 



air: 



0.29m^ 



(A.l) 



ai = 0.286 -f 0.626^6 

a2 = 0.1786 - 6.308wfc + 174.9cjg-1168cjg (A.2) 



<72 = 0.56[1 + 21 Awl^Y^ ■ (A.6) 

The relative shift of the first acoustic peak is zero, 
Sifii — 0, and the relative shifts of the second and third 
peaks are given by 



S(p2 — Co - cir* - C2r^ -t- 0.05(ns — 1) 



(A.7) 



with 



0.1+ 0.213-0.123^2 



X exp |- (52 - 63.6f^t) ujb^ , 

ci = 0.015 -I- 0.063 exp (-3500wg) , 
C2 = 6 X 10"^ + 0.137(cjf,-0.07)2, 

C3 = 0.8 + 2.3nt + (70 - 126J7f^) ojb , (A.8) 

and 

Sip3 = 10 - dirf^ + 0.08(n, - 1) , (A.9) 

with 

di = 9.97+ (3.3 - LUb, 

rfa = 0.0016 - 0.006717^ + (0-196 - 0.22^^^ ojb 

+ (2.25 + 2.77Cll^ X lO~^uj-\ (A.IO) 

The relative shift of the first trough is given by 



Sy:'3/2 ^bo + birl^^ eMhr,) + 0.158(71^ - 1) (A.U) 
with 
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The estimate for the accuracy of the fitting formulae 
_ is given at la by [20| : 

bo = -0.086 - O.Orgr^Jj* - (2.22- 18. If7^")cjf, 

- (140 + 403^'^) Lof , 

^ ^' . A,^ = 0.0031, 

b, = 0.39 - 0.98^^^ - (iS.l - 29.2ny ^^^^ ^ OO44 

+UOujI , AS(j)y2 = 0.0039 . (A. 13) 

62 = -0.57- 3.8 cxp(-2365cj^) . (A.12) 
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